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Abstract 

The evolution of matter in the expanding FRW uni- 
verse during the time interval between the end of in- 
flation and the beginning of the radiation-dominated 
era is studied. A constraint between the global ge- 
ometry and total amount of matter in the universe 
as a whole, which is valid during the phase of an in- 
tensive transfer of energy to the matter degrees of 
freedom, is introduced. The matter is considered as 
a perfect fluid with two components between which 
there is energy exchange. The analytical solutions of 
the Einstein equations are found. The limiting cases 
of the the Hubble expansion rate and the total en- 
ergy density, which correspond to matter production, 
pressure-free and radiation-dominated phases are in- 
vestigated. The transition to the inflationary phase 
and a unidirectional evolution of matter in the uni- 
verse at all phases are discussed. 



1 Introduction 

The standard cosmological model is based on the hot 
Big Bang model and the assumption about the in- 
flationary expansion of the very early universe, when 
the scale factor grows quasi-exponentially, while the 
Hubble expansion rate remains almost constant (e.g., 
Refs. [I] [51 121 H] ) . During inflation the universe is in 
the vacuum-like state which is usually associated with 
a scalar field called the inflaton. After inflation the 
energy density of the primordial matter (except the 



inflaton) which filled the universe before this stage 
becomes negligibly small. In order to explain the 
presence of conventional matter in the universe af- 
ter inflation, the decay of the vacuum-like state into 
'normal' particles is postulated. 

The universe becomes hot as a result of interac- 
tion between particles and transits into the radiation- 
dominated phase. In the process of energy transfer 
from the inflaton to radiation (called reheating) the 
equation of state of matter changes. A change from 
a vacuum-like equation of state to the equation of 
relativistic matter might be gradual and an interme- 
diate stage between these two known phases may be 
modeled. 

In the present article, the evolution of the equa- 
tion of state of the matter in the universe during the 
time interval between the end of inflation and the be- 
ginning of the radiation-dominated era is considered. 
Without rendering concrete mechanisms of decay of 
vacuum-like state into the conventional matter, we 
assume that the global geometry and total amount of 
matter in the universe as a whole satisfy a constraint, 
which is valid during some time interval, before radia- 
tion domination. This constraint is equivalent to the 
law of the conservation of total energy of the universe 
which remains equal to zero due to the gravitational 
mass effect, whereas the energy attributed to the par- 
ticles of conventional matter increases with expansion 
of the universe [H [5] . In this case, at all stages of 
evolution the universe is described by the Einstein 
equations with addition of appropriate equations of 
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state. 

The paper is organized as follows. In Section 2 
the basic equations which describe the homogeneous, 
isotropic and spatially flat universe are given. The 
equations of state of matter for the different phases of 
reheating are justified. In Section 3 a two-component 
perfect fluid model is introduced. The analytical so- 
lution of the non-linear equation for the Hubble ex- 
pansion rate is obtained. The expressions for the de- 
celeration parameter and the total energy density as 
functions of time are deduced. The limiting cases 
of the solutions which correspond to pressure-free 
and relativistic matter are considered. The Whitrow- 
Randall's relation [B] is rederived. In Section 4 the 
transition to the inflationary phase is discussed. The 
mechanical analogy which explains a unidirectional 
evolution of matter in the universe at the phases un- 
der consideration is given. 



2 Equation of state parameter 

Let us consider the homogeneous, isotropic and spa- 
tially flat universe in the early epoch, when its dy- 
namics can be described by the equations of the 
Fricdmann-Robertson- Walker (FRW) cosmology, 



H 2 = 



R 
R 



8irG 



-p, 



p + 3H(p+p) = 0, 
V = w(t) p, 



(1) 



(2) 
(3) 



where R(t) is the cosmic scale factor, pit) is the en- 
ergy density of the matter which has a form of the 
homogeneous perfect fluid, p{t) is its pressure, w(t) is 
the equation of state parameter, G is the Newtonian 
gravitational constant, an overdot denotes d/dt, t is 
the proper time (units c = 1 arc used). 

Let there exist an interval of time after inflation 
A£ c i = t c — ti, where U denotes the time at which 
inflation ends and t c stands for the time at which an 
intensive transfer of energy to the matter degrees of 
freedom ends. We assume that during this interval 



the matter is produced, so that the following condi- 
tion is fulfilled, at least in good approximation, 



M 2 

M - G— = 0, 



(4) 



where M — ^nR?p is total mass-energy of matter 
(the sum of masses of particles of conventional mat- 
ter) in the equivalent flat-space volume taken with- 
out account of gravitational interaction between par- 
ticles. The equation Q can be interpreted as the 
law of the conservation of zero total mass-energy of 
the universe during its expansion with matter pro- 
duction 5 . According to Eq. (QJ, during the time 
interval At C i the following relation R = GM is valid. 
It means that the energy density 
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GA-kR 2 



(5) 



decreases linearly with increasing surface area AirR 2 . 

From Eqs. ([2) and one can find the equation 
of state | | 

P=-gP, ™(t) = -~. (6) 

After the end of this phase, the mass M remains 
constant on the time interval Ai rc = t r — t c , where t r 
denotes the beginning of the subsequent relativistic 
matter dominant era. The equation of state on the 
time interval At rr takes the form 



p = 0, wit) = 0. 



(7) 



In the relativistic matter dominant era, the mass 
attributed to relativistic matter reduces as the uni- 
verse expands, M ~ R^ 1 , due to the cosmic redshift. 
For times t > t r , the equations of state has a form 



1 . . 1 

P=gP, w{t)=-. 



(8) 



We will study the model of evolution of matter in 
the early universe, where the equation of state pa- 
rameter wit) changes with time from — i, passing 
through the point w = 0, to 4 taking all intermedi- 



ate values, 
of time t, 



Substituting a continuous function wit) 



wit) 



tanh 



t - tr 



(9) 
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for the equation of state parameter on the time inter- 
val [ii,t r [, and choosing properly a point to on this 
interval, one can reproduce the required values of 
w ([5])-©. Since t r ^> ti (in standard cosmological 
model, the value fj ~ 10 -35 s is acceptable, whereas 



the time t r is often evaluated as t r 



10 



-30 



responding to temperatures not exceeding 10 12 GeV 
PQ), the good estimation for t may be to <t r . The 
value 1 /r determines the mean rate of change of the 
equation of state parameter w(t). Such a variation 
of the equation of state parameter can be achieved in 
a system, where the matter consists of a few compo- 
nents between which occurs the energy transfer for 
some typical time r. 

3 Two- component fluid 

Let us consider a two-component perfect fluid with 
the energy density and pressure 

P = Pq+Pd, p=p q +Pd- (10) 

These components satisfy the equations 

p q +3H(p q +p q ) = Q, p d +'SH(p d +Pd) = -Q, (11) 

which represent the energy conservation law ([2]) 
rewritten for components, Q is the interaction term. 

The components of the perfect fluid are imitated 
by scalar fields <fr q (t) and <j>d(t) with potentials V q {<fi q ) 
and V d {(/)d), 
1 

Pa = V 

The models of such a type which include a coupling 
between the matter components were considered in 
the literature, in particular, within the context of in- 
flation and reheating and the coincidence problem of 
dark energy and matter in the accelerating universe 
(see, e.g., Refs. [71 [51 IS] and references therein). The 
form of the interaction term Q may be derived from 
different physical arguments or obtained as a solution 
of some dynamical equation, which describes the re- 
quired properties of the matter fields (fi a . 

Let us assume that the field (fid forms the pressure- 
free matter component (dust), 



2 1*2 

« + V ai p a = -<fi a - V a , a = {q, d}. (12) 



1 • 2 

-ifi d = Vd, Pd = 2Va, 



Pd 



0. 



(13) 



Concerning the field <fi qi we suppose that it is de- 
scribed by the vacuum-type equation of state (as for 
the inflaton) at times t <C t , 



Pq 



(14) 



From Eq. (|12p. it follows that at this stage the kinetic 
energy of the field <fi q can be neglected and the total 
energy is determined by its potential term, 



0, Pq ~ V q . 



(15) 



For times t 
form 

It means that 



to, the equation of state takes the 

Pq * 0. (16) 



n^q 2 



V q , Pq ~2V q . (17) 

Then, for the times t 3> to, the field <fi q describes the 
matter component with the energy density which is 
almost equal to its kinetic energy, 



Pq 



V a ~ 0. 



(18) 



This phase corresponds to the reheating of the 
pressure-free matter and provides the passage to rel- 
ativistic matter domination. The field <fi q here has a 
form of the stiff Zel'dovich matter, 

Pq^pq at t > t . (19) 

The continuous transition from Eq. (fT4")l to ([TB"]) . and 
then from ([TO]) to ([T9]) can be achieved if the following 
condition is imposed on the field <fi q 



\<t>q 



e -2(t-i )/r = y 



(20) 



where r < ^to. 

Taking into account Eqs. JTDJ), (JT3J), dTSJ), 

PH). from Eq. (J3J we find 



w(t) 



e 2(t-t )/r _ i 



(21) 



e 2(t-t )/r + 1 + 2V d /V q ' 

This relation passes into Eq. ©, if one introduces 
the following additional condition on Vd, 



V d = p q = V q 



i(t-to)/r 1 



(22) 
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Then from Eqs. (JD) and ([13]), we get 



P = 3p q , P = Pq, w 



Pq_ 
3Pq' 



(23) 



In this case, the interaction term Q — 2Hp q and the 
set of equations (ITT1) reduces to one equation 



(24) 



From Eqs. ©, ©, and HU, it follows the non-linear 
equation for the Hubble expansion rate, 



H + ~ <! 3 + tanh 



t-tr 



\ I t/2 



ff^ = 0. (25) 



The general solution of this equation is 

2 



H(t) = 



D(ty 



where we denote 



D(t) = Ct + 3t + rln cosh 



t-t 



(26) 



(27) 



C is a constant of integration. 

The deceleration parameter, q = — 1 — H/H 2 , is 
equal to 



?(*) = ^ M 



tanh 



t - if 



(28) 



The deceleration parameter changes from the value 
q = for the equation of state ©, through the point 
q = \ for Eq. ©, to g = 1 for Eq. ©. Thus, 
in the model under consideration, the expansion of 
the universe is decelerating on the whole time inter- 
val from the end of inflation to the beginning of the 
radiation-dominated era. 
The total energy density is 



Pit) 



2nGD(ty 



(29) 



The limiting cases of the solutions (1261) and (|29| 
reproduce the well-known expressions for the Hub- 
ble expansion rate and the energy density. Setting 



C = 0, near the point t = to we find for pressure- free 
matter [10] 

(30) 

Choosing the constant C ~ 1 + ^- In 2, for t 3> to > 2r 
we obtain the relations for the relativistic matter 



ff ^5' ^32^ 



(31) 

For times t -C to and to > 2r, from Eq. ([27]) it follows 
D(t) = 2t+(C+l)t -rln2. (32) 
Setting C~ — 1 + ^ ln2, these expressions reduce to 

1 3 

(33) 



^4 



The equation for p(t) has a form of Whitrow- 
Randall's relation [6]. 

4 Discussion 

The equations ([2"5)) , ([2T]) . and ([2^]) demonstrate how 
the Hubble expansion rate and the energy density 
change with time from the inflationary phase of the 
universe's expansion, through the subsequent eras of 
an intensive energy transfer and pressure-free mat- 
ter, to the beginning of the radiation domination. By 
choosing the constant of integration C, the solutions 
([26]) and ([29]) are reduced to known 'standard' ex- 
pressions ([3"0"]) . ([5T]) . and ([3"3"]) . An interesting feature 
of the solution ([29l) is that at the point t = it is 
finite, 



P(0) 



27rG[(C* + l)< -Tln2] 2 ' 



(34) 



Thus, the two-component system does not have an 
initial cosmological singularity. 

The equations ([26]) . ([29]) . and ([32]) can be continued 
into the region of extremely small values of time, t <S 
1 1 (C + 1) to — t In 2 1 , where the Hubble expansion rate 
slightly changes with time, so that in the inflationary 

phase H{ti) ~ H(0) = J ^p/?(0) and the expansion 
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of the universe will be exponential in time, R(t) ~ 
exp{H(ti)t}. 

The expression for the energy density p(U) in the 
inflationary phase can be reduced to the 'standard' 
form. Setting G = Alp 2 [3J, where Mp is the Planck 
mass, t ~ Mp , and choosing C = — 1, from Eq. 

with a good accuracy we get H(ti) ~ yJ^Mp, 

p(U) ~ Mp 4 (cf. P). 

In the two-component model (fTU)) with the equa- 
tion of state ([3]) with the parameter , the evolution 
of the universe goes in one direction, from small times 
t <C to to large values t 3> to- 

The following mechanical analogy allows to under- 
stand the reason of the origin of this 'arrow of time'. 
The function ^ can be considered as the kink solu- 
tion of the equation 

ltf + [-U(w)]=0, U = ^^vr i -^j , ( 35 ) 

which describes the motion of the analogue particle 
with zero energy in the potential [— U(w)} (cf., e.g., 
Ref. [II])- This potential has two maxima at the 
points w = ±jj and a local minimum at w = 0. 
The analogue particle moves along the 'trajectory' (JSJ 
from the value w = — -| in the distant past (t = — oo) 
to the value w = | reached at t = oo. At the moment 
t = to i the particle passes through the minimum of 
the potential at w = 0. Leaving the point w = — h, 
the analogue particle can only approach the point 
w — h at t — > oo, where its velocity and acceleration 
vanish. It cannot return back to w = — h ■ 

We note that Eq. (|35|) has another solution in the 
form of the antikink which is equal to the function 
<j9j) with an inverse sign. This case corresponds to 
the model in which the relativistic matter at t = — oo 
transforms into the pressure-free matter and then 
into a gas of low- velocity cosmic strings at t = oo. 
It was studied in Ref. [12], where it was shown that 
the equation of state of matter can change with the 
expansion of the universe due to energy transfer be- 
tween the matter components (scalar fields) allowing 
to reproduce the evolution of matter in the universe 
with non-zero cosmological constant. The description 
on equal footing of the universe over the total time in- 
terval from inflation through reheating to subsequent 



cooling and transition to the pressure-free matter us- 
ing the kink and antikink solutions of Eq. (1351) may 
indicate about their important role in the represen- 
tation of the evolution of matter. 
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